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We investigate the metal-insulator transition in the half-filled Hubbard model on a two-dimensional 
triangular lattice using both the Kotliar-Ruckenstein slave-boson technique, and exact numerical 
diagonalization of finite clusters. Contrary to the case of the square lattice, where the perfect 
nesting of the Fermi surface leads to a metal-insulator transition at arbitrarily small values of U, 
always accompanied by antiferromagnetic ordering, on the triangular lattice, due to the lack of 
perfect nesting, the transition takes place at a finite value of U, and frustration induces a non- 
trivial competition among different magnetic phases. Indeed, within the mean-field approximation 
in the slave-boson approach, as the interaction grows the paramagnetic metal turns into a metallic 
phase with incommensurate spiral ordering. Increasing further the interaction, a linear spin-density- 
wave is stabilized, and finally for strong coupling the latter phase undergoes a first-order transition 
towards an antiferromagnetic insulator. No trace of the intermediate phases is instead seen in 
the exact diagonalization results, indicating a transition between a paramagnetic metal and an 
antiferromagnetic insulator. 
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The Mott metal-insulator transition (MIT), i.e. the 
transition from a metallic to an insulating phase driven 
by the electronic correlation jl],^] , is one of the most rel- 
evant issues in condensed matter theory. In the last few 
years it has been also the object of an intensive study, 
due to many experimental evidences of Mott insulators 
ranging from the parent compounds of the superconduct- 
ing cuprates Q to the alkali fullerides of the type A4C60 

i- 

The simplest model in which the competition between 
the delocalizing effect of the kinetic energy and the elec- 
tronic correlation can give rise to a MIT is the Hubbard 
model 



n 
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<ij> 



(4. 



h.c 



(1) 



where c] a {c^ a ) creates (destroys) an electron with spin a 

on the site i and Ui, ff = c\ a a^ is the number operator; 
t is the hopping amplitude, U is the Hubbard on-site 
repulsion, fi is the chemical potential. The hopping is 
restricted to nearest-neighbors and the indices i,j label 
the points and Vj of a d-dimensional lattice. 

At half-filling (i.e. for a number of electrons equal to 
the number of sites), this model is known to undergo a 
MIT by increasing the interaction strength U. On the 
d-dimensional cubic lattice the perfect nesting property 
of the Fermi surface makes the model unstable towards 



antifcrromagnetism as soon as a non-zero U is turned on, 
driving the system to the insulating state. In this paper 
we focus on the triangular lattice as a prototype for a 
model where the perfect nesting is absent for the uncor- 
rected metal ||. Since in the U/t — > 00 (Heisenberg) 
limit the model is likely to display a Neel ordered (insu- 
lating) ground state (GS) Hf7|, the MIT is expected to 
occur for finite U . 

Besides its theoretical relevance, our analysis has also 
an experimental counterpart. In fact, the adlayer struc- 
tures on semiconductor surfaces, such as SiC(OOOl) || or 
K/Si(lll) : B Q, have recently turned out to be an al- 
most ideal environment for the study of Mott insulators 
|lfjf and are characterized by a y/3 x \/3 arrangement 
of the dangling-bond surface orbitals, which are likely to 
be well described by bidimensional strongly correlated 
Hamiltonians jnj on the triangular lattice. 

Hartree-Fock (HF) calculations, performed by Krish- 
namurthy and co-workers [ p~T|Jl2| , produce a rather rich 
phase diagram: for small U the system is a paramag- 
netic metal (PM), which turns to a metal with incom- 
mensurate spiral spin-density-wave (Spiral Metal, SM) at 
U = U c i — 3.97i. Two successive first order transitions 
occur further increasing the coupling: at U = U C 2 = 4.45< 
a semi-metallic linear spin-density wave (LSDW) is sta- 
bilized, and a first order MIT to an antiferromagnetic in- 
sulator (AFMI) occurs at U = U C 3 = 5.27t. In the same 
work it has also been argued that at finite temperature 
the model should present a Mott transition between a 
paramagnetic metal and a paramagnetic insulator. 
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However, the above transitions only occur at relatively 
large U/t and the HF approximation is unreliable in the 
intermediate- and strong-coupling regime. Therefore we 
adopt the more appropriate slave-boson (SB) approach 
JT3| , fl4[ as an interpolating scheme between the U/t = 
and the U/t — > oo regimes. To allow for the presence 
of incommensurate spiral spin ordering, we introduce the 
spin-rotational invariant formulation Jl5| of the Kotliar- 
Ruckenstein SB approach iQ. The reader can find fur- 
ther details in Ref. Jl5[ |. 

We introduce on each site a set of four SB operators 
ej, Sj,-f) and di to label empty (e), singly (s), and dou- 
bly (d) occupied sites, respectively. The spin projection 
<T =|, I is measured with respect to a local quantization 
axis, which is allowed to vary from site to site. The re- 
sulting SB Hamiltonian is 

<ij> 

-vJ2 d L d ^ + u Y, d l d * 

+ E ^eje, + d\di + 4,, s i* ~ 1 j 

where Ci i? ,cJ f are the pseudofermion operators, the La- 
grange multipliers Xi and Ai jS enforce on each site the 
completeness constraint and the correct fermion counting 
respectively, the operator IZi rotates the local reference 
frame back to the laboratory frame, and the operator 



(2) 



1 - d\di - s\ s s is V 1 - 4 e i ~ 4 



reconstructs the hopping amplitude in the enlarged Fock 
space, and yields the correct U — * limiting behavior 
in the mean- field approximation |lj,|l5|. When the an- 
gle between two local quantization axes depends only on 
their relative position, up to a global phase factor one can 
assume 'R-l'R-j — exp[iQ • (Rj — R J -)r 3/ /2], where r y is the 
Pauli matrix and Q is the (incommensurate) modulating 
wavevector p"5[ . In such a case a mean- field description 
with real site-independent expectation values for the SB 
operators 



(e 



(t)\ 



eo; 



„(t)\ 



so,. 



R (t) > 



(3) 



and for the Lagrange multipliers 



(A l )=A ; (A i ,,)=A ,„ (4) 

is possible. Eqs. (^) and (^) refer to the case in which 
the translational symmetry is not broken and the expec- 
tation values of the bosons and of the Lagrange multi- 
pliers do not depend on the site. We have also studied 



configurations with broken translational symmetry. In 
particular we considered solutions in which in the bosons 
have different values on each of the three sublattices, and 
analogous to LSDW found in Hartree-Fock jll] . The lat- 
ter solutions can be found considering a four-sites unit 
cell. A similar SB calculation has been performed in Ref. 
]l6[ , where, however, the generalization of the SM phase 
found in HF was never recovered as an energy minimum. 
In the case of spiral spin ordering, the Hamiltonian (|^) 
can be analytically diagonalized by adopting the Bloch 
representation, and performing a unitary transformation 
with respect to spin indices, yielding 



Eh 



1 



k,± 



± 



t{zf 



o,T 



A 0:T - A 



o.l 



Af2 7 2 7 1 rp2 

%i 2 ,T z o,| J o 



where T e = - J2i cos (Q ■ V 2 ) cos(k • I), T a = - sin(Q • 
1/2) sin(k-l), and 1 = (1,0), (1/2, ±V§/2) are the nearest- 
neighbor displacements. The self-consistency equations 
are obtained by minimizing the free energy 



T = T a -T lo § 



1 



-E Ka /T 



whereto = N[U(^+X (eo+'^+ s o,'t +s o,l~ 1 )~ A <i,'t(^+ 
Sq l) — A 0i |(c?§ + Sq |) + fin], N is the number of sites, 
and n is the electron density per site, and read 



9^0 \ - 
dX ^ 



dX 



-/(£*,«) = 0, 



(5) 



where f(E) = [e E ' T + is the Fermi function and X 
represents generically one of the parameters (^), (Q) and 
the two components of the pitch vector Q. The chemical 
potential fi is fixed by the condition 



^ f{E Ka ) = nN. 



k,a=± 

In this paper we assume henceforth n = 1 (half-filling). 

The self-consistency equations (|^) yield the same so- 
lutions found in HF, namely a paramagnetic metal, a 
metal with incommensurate spiral ordering, a linear spin- 
density-wave and an antiferromagnetic insulator. As in 
HF, the PM-SM transition is continuous, and the other 
two transitions are of first order, but all of them occur 
at larger coupling values, U c \ = 6.68i, U c i — 6.84f, and 
£/ C 3 = 7.68t. The energy curves corresponding to the 
above phases are reported in Fig. [I]. Our results agree 
with Ref. @ as far as the PM, AFMI, and LSDW phases 
are concerned, but we also find a region of stability for the 
SM phase, which was not detected in Ref. . These au- 
thors were indeed looking for spiral phases starting from 
the strong-coupling side, and following them to weaker 
coupling. On the other hand, our analysis shows that a 
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spiral phase develops continuously from the PM at in- 
termediate coupling and it ends in a critical point soon 
after the level-crossing with the AFMI (see the inset in 
Fig. [j]), and does not exist at strong-coupling. Therefore, 
our SM phase is the generalization of the corresponding 
phase found within HF |L2|, and it is unrelated to the 



high-energy SM phases of Ref. |16|. However, the region 
of existence of the SM is narrower within SB as compared 
to HF, and the magnetization m = \{n-\ — njj is always 
less than 0.1, a really small value with respect to the 
HF value (< 0.4). Therefore the jump of the magneti- 
zation at the SM-AFMI transition is substantially larger 
than in the HF approximation. We point out that, con- 
trary to nesting models, where the presence of free parti- 
cles (doping) is a necessary condition for spiral ordering 
^5| , p7|Jl8| , here the spiral phase exists at half-filling, as 
previously shown in Ref. |Tl| l , within the HF approxima- 
tion. Despite the overall qualitative agreement between 
the HF and the SB phase diagrams, the main outcome 
of the comparison between them is that the stability of 
the SM phase is strongly reduced. Furthermore, the SM 
is hardly distinguishable from the PM in its whole region 
of stability. It is reasonable to expect that the inclusion 
of quantum fluctuations washes out these phases leav- 
ing the way open for a transition between a PM and the 
AFMI. 
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FIG. 1. U dependence of the ground-state energy per 
site. SB results: PM (solid line), SM (dotted line), LSDW 
(dot-dashed line), and AFMI (dashed line). Open dots are 
the exact diagonalization results for the N = 12 cluster. 

Despite the strong frustration of the antifcrromagnctic 
(AFM) order on the triangular lattice §§], both the HF 
and SB approaches indicate no paramagnetic Mott insu- 
lating phase in the zero-temperature phase diagram of 
the half-filled Hubbard model. 

In particular, within the SB approach, we can indicate 
how far the system is from the Brinkman-Rice transi- 
tion Q to a paramagnetic Mott insulator. In fact, if 
the possibility for magnetic ordering is neglected, the 
paramagnetic metallic phase undergoes a Brinkman-Rice 



transition with vanishing double occupancy and effective 
hopping amplitude, for a critical value (at T — 0) of the 
Hubbard interaction Ubr = 32£iV~ 1 ^ k £k@(2tek + m) 
(c* = -T e (Q - 0)) |Q, i.e. U BR /t ~ 15.8 on the 
triangular lattice. As we see this value is much larger 
than U c i/t, U c ijt and U c s/t found above. The system is 
therefore not even close to the Brinkman-Rice transition 
when the MIT occurs. 
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FIG. 2. Overlap of the finite-!/ GS with the U = (empty 
dots) and the U = WOt (full dots) GS's, for TV = 12. In the 
inset the GS energy per site in the k = (0, 0) (solid line) and 
k = (2n/3, 0) (dashed line) subspaces is plotted versus U. 

In order to understand to which extent the picture we 
found within the mean-field SB theory survives in an ex- 
act treatment of the model, we performed exact diag- 
onalization of small clusters by means of the standard 
Lanczos algorithm. The largest lattice compatible with 
all the symmetries of the model that can be handled with 
exact diagonalization is a N — 12 site cluster ||. We al- 
ways used twisted boundary conditions with a suitable 
phase such that the half-filled system is in a closed-shell 
configuration. This is important in order to perform a 
reasonable investigation of the conduction properties of 
the finite-size system. It turns out that the boundary 
conditions that minimize the energy in a closed-shell con- 
figuration for U — leave the system in a closed-shell 
configuration at all U. The energy is shown as a function 
of U in Fig. [l[ The overall agreement with the mean-field 
SB results is good, the largest deviations (~ 20%) being, 
as expected, at intermediate coupling (U/t ~ 7). 

To check the occurrence of a discontinuous phase tran- 
sition we evaluated the overlap between the GS wave 
function and the two limiting cases of U = 0, and for large 
U (namely, U — lOOt). As shown in Fig. ^ on the large- 
U side of the diagram the GS has a large overlap to the 
AFM strong-coupling state and a vanishing overlap with 
the non-interacting metallic one. On the metallic side the 
overlap with the non-interacting state is always finite, but 
it is a decreasing function of U; in this regime the GS has 
anyway a vanishing overlap to the AFM state. We have 
therefore a clear evidence for a strongly correlated metal 
with a decreasing coherent part. In particular the sharp 
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change of the GS wave function at Umit — 12.07t is due 
to a level-crossing occurring between a metallic and an 
antiferromagnetic solutions, as it is shown in the inset 
of Fig. ||. These results, however, do not rule out the 
possibility of a continuous transition within the metallic 
phase, i.e. the PM-SM transition found with SB. 
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FIG. 3. Spin structure factor 5(q): U=2t (dotted 
line), U=8t (dashed line), U=11.5t (dot-short-dashed line), 
U=12.5t (dot-long-dashed), U=40t (solid line). 

In Fig. H we show the spin structure factor, 5(q) = 
J2i j SiS* exp[q • (r$ — Tj)]/N, for different values of U. 
The results do not suggest any intermediate state be- 
tween a metallic state without magnetic order and the 
AFM insulator, as 5(q) abruptly changes from a struc- 
tureless behavior to an AFM pattern peaked at the clas- 
sical ordering wavevector, i.e. Q = (47r/3, 0). Although 
we suspect that the intermediate phases are an artifact 
of the mean-field approach, the weakness and the strong 
size-dependence of the spiral phases suggested by the SB 
results, may make them unaccessible on our 12-site lat- 
tice. 




FIG. 4. Exact calculation of the Drude weight as a function 
of U for the N = 12 cluster. 

Using the Lanczos algorithm we have also calculated 
the finite-frequency optical conductivity <j(u>) and the 
Drude weight, measuring the electronic mobility. The 



real part of the xx component of the conductivity tensor 
for a tight-binding model at zero temperature may be 
expressed in terms of the Kubo formula |nj] 

o»*M = D xx 6(cj) + 3<0| 4 — A- -J x |0), (6) 

uj — n + Lq — id 

where J x = J2i,a,i ^( c l,<r c «'+ 1 >< r - h - c -) is tlie ^-component 
of the current operator. The coefficient of the zero- 
frequency delta function contribution D xx , the Drude 
weight, is given by the f-sumrule Q 



D, 



2 {n * } ^ E n -E ■ 



(7) 



where U x = Y,i, cr ,i l l{ c i,<T C i+^ + h - c -)> and I0n) is tlic 
eigenfunction of Ti. with eigenvalue E n . 

The latter quantity, which is reported in Fig. |] is a 
direct measure of the metallic character of the state, and 
the MIT is signaled by the vanishing of D xx (^0|. For 
a finite system, D xx does not vanish for any value of U, 
but an abrupt change takes place at the level-crossing 
point. For U < Umit, Dxx is a decreasing function of 
the interaction, which resembles the overlap in Fig. |^. 
An abrupt change takes place at the level-crossing point 
and for U > Umit it becomes negative, a common phe- 
nomenon in the insulating phase of a small-size system 

All the results of exact diagonalization point towards 
the same direction: the metal- AFMI level-crossing found 
within the SB mean-field approach is shifted to larger 
values of U. The metallic solution exhibits a continuous 
loss of metallicity with increasing U. The Drude weight 
is finite up to the MIT on the 12-site lattice although 
it is quite small (4% of the non-interacting value). We 
remark that, due to finite-size effects, we cannot exclude 
the possibility that D xx vanishes before the transition to 
the AFMI is reached. In such a case, there would be 
a region of parameters in which the paramagnetic insu- 
lator exist, though the SB results point in the opposite 
direction. 

In conclusion, using the slave boson technique and 
the exact diagonalization, we have investigated the zero- 
temperature phase diagram of the half-filled Hubbard 
model on a two dimensional triangular lattice. The 
mean-field SB approach displays a rich phase diagram 
which qualitatively resembles the one from HF calcula- 
tions, but, on the other hand, drastically reduces the 
stability of the spiral metal and of the linear spin- 
density-wave states. Namely, the weak-coupling para- 
magnetic metal continuously evolves into a spiral metal 
at U = U c i = 6.68i, which crosses the linearly polar- 
ized spin-density-wave ground-state at U = U c i = 6.84£. 
The latter phase undergoes a further first-order tran- 
sition towards an antiferromagnetic insulator at U = 
U C 3 = 7.68t. All these transitions occur for coupling con- 
stants substantially smaller than the critical value for the 
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Brinkman-Rice transition to a paramagnetic insulator 
(Ubr = 15. 8t). The exact-diagonalization results present 
a first order transition between the paramagnetic metal 
and the antiferromagnetic insulator at Umit = 12.07i, 
without intermediate "exotic" phases. 
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